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<N ! Abstract 



Numerical simulations indicate that the smooth, photoionized intergalactic medium 
(IGM) responsible for the low column density Lyo forest follows a well defined temperature- 
density relation. We show that such an equation of state results in a cutoff in the distribu- 
tion of line widths (6-parameters) as a function of column density (N) for the low column 
density absorption lines. This explains the existence of the lower envelope which is clearly 
seen in scatter plots of the 6(iV)-distribution in observed QSO spectra. The intercept and 
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, slope of this cutoff can be used to measure the equation of state of the IGM. Measuring 



the evolution of the equation of state with redshift will allow us to put tight constraints 



O . on the reionization history of the universe. 
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22 1 1 Introduction 
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Neutral hydrogen in the intergalactic medium (IGM) gives rise to a forest of Lya absorption lines 
blueward of the Lya emission line in quasar spectra. Hydrodynamic simulations of structure 
formation in a universe dominated by cold dark matter and including an ionizing background, 
show that the low column density (N ^ 10 14 ' 5 cm -2 ) Lya forest at redshift z ^ 2 is produced by 
a smoothly varying IGM. For the low-density gas responsible for the Lya forest, shock heating 
is not important and the gas follows a tight temperature-density relation. The competition 
between photoionization heating and adiabatic cooling results in a power-law 'equation of state' 
T = T (p/p) 7_1 , which depends on cosmology and reionization history ||. 

The smoothly varying IGM gives rise to a fluctuating optical depth in redshift space. Many 
of the optical depth maxima can be fitted quite accurately with Voigt profiles. The distribution 
of line widths depends on the initial power spectrum, the peculiar velocity gradients around the 
density peaks and on the temperature of the IGM. However, there is a lower limit to how narrow 
the absorption lines can be. Indeed, the optical depth will be smoothed on a scale determined 
by three processes: thermal broadening, baryon (Jeans) smoothing and possibly instrumental, 
or in the case of simulations, numerical resolution. The first two depend on the thermal state of 
the gas. While for high-resolution observations (echelle spectroscopy) the effective smoothing 
scale is not determined by the instrumental resolution, numerical resolution has in fact been 
the limiting factor in many simulations (see || for a discussion). 

Scatter plots of the 6(A)-distribution have been published for many observed QSO spectra 
[[[], [6|, ||, f|. These plots show a clear cutoff at low 6-parameters, which increases slightly with 
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Figure 1: The &(iV)-distribution for 800 random lines from models with log To = 4.20 (a) and 
log To = 3.83 (b) respectively, at z = 3. The position of each line is indicated by a cross. Errors 
are not displayed. The dashed line is the cutoff for the lines plotted in panel (a) over the range 
10 12 - 5 cm" 2 < N < 10 14 - 5 cm" 2 . 

column density. However, this cutoff is not absolute, there are some narrow lines, especially 
at low column densities. Lu et al. || and Kirkman & Tytler || use Monte Carlo simulations 
to show that many of these lines are caused by line blending and noise in the data. Some 
contamination from unidentified metal lines is also expected. A lower envelope which increases 
with column density has also been seen in numerical simulations ||. 

In this contribution we shall demonstrate that the cutoff in the 5(iV)-distribution is deter- 
mined by the equation of state of the low-density gas and can therefore be used to measure 
the equation of state of the IGM. This work will be more fully described and discussed in a 
forthcoming publication |7J. 

2 Results 

In Fig. ||a we plot the 6(iV)-distribution for 800 random absorption lines taken from multi- 
component Voigt profile fits of spectra at redshift z = 3, generated from one of our simulated 
models. A cutoff at low 6-values, which increases with column density, can clearly be seen. As 
in the observations, there are some very narrow lines, which occur in blends. These unphysically 
narrow lines make determining the cutoff in an objective manner nontrivial. We developed an 
iterative procedure for fitting a power-law, b = b^N/No^^ 1 , to the 6(iV)-cutoff over a certain 
column density range (10 12 5 cm -2 < N < 10 14 ' 5 cm -2 at z — 3) which is insensitive to these 
narrow lines. 

The 6(iV)-distribution for a colder (log To = 3.83 vs. log To = 4.20) model is plotted in 
Fig. ^b. Clearly, the distribution cuts off at lower 6-values. Let us assume that the absence of 
lines with low b- values is due to the fact that there is a minimum line width set by the thermal 
state of the gas through the thermal broadening and/or baryon smoothing scales. Since the 
temperature of the low-density gas responsible for the Lya forest increases with density, we 
expect the minimum 6-value to increase with column density, provided that the column density 
correlates with the density of the absorber. 

To see whether this picture is correct, we need to investigate the relation between the Voigt 
profile parameters iV and b, and the density and temperature of the absorbing gas respectively. 
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Figure 2: The gas density at the line centres as a function of neutral hydrogen column density. Data 
points correspond to the lines plotted in Fig. [l]a. 
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Figure 3: The temperature at the line centres as a function of 6-parameter. The left panel contains 
all lines plotted in Fig. [l|a. In the right panel only those lines are plotted which have 6-parameters 
within one mean absolute deviation of the power-law fit to the 6(iV)-cutoff plotted as the dashed line 
in Fig. ||a. The dashed line corresponds to the thermal width, b = (2A;sT/m p ) 1 / 2 . 

From Fig. 0, which shows the gas density as a function of column density for the absorption 
lines plotted in Fig. [l]a, it can be seen that these two quantities are tightly correlated. 

The temperature is plotted against the 6-parameter in Fig. |3]a. The result is a scatter plot 
with no apparent correlation. This is not surprising since many absorbers will be intrinsically 
broader than the local thermal broadening scale. In order to test whether the cutoff in the 
6(iV)-distribution is a consequence of the existence of a minimum line width set by the thermal 
state of the gas, we need to look for a correlation between the temperature and 6-parameters 
of the lines near the cutoff. Fig. |]b shows that these lines do indeed display a tight correlation. 
The dashed line corresponds to the thermal width, b = (2kgT ' /m p ) 1 / 2 , where m p is the mass of 
a proton and ks is the Boltzmann constant. Lines corresponding to density peaks whose width 
in velocity space is much smaller than the thermal broadening width, have Voigt profiles with 
this 6-parameter. 

Figs. H] and |]b suggest that the cutoff in the 6(iV)-distribution should be strongly correlated 
with the equation of state of the absorbing gas. The objective is to establish the relations 
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Figure 4: The intercept of the 6(iV)-cutoff as a function of the temperature at mean density (left 
panel) and the slope as a function of the index of the power-law equation of state (right panel), for 
redshift z = 3. The error bars enclose 68% confidence intervals around the medians, as determined 
from 500 sets of 300 random lines. The dashed lines are the maximum likelihood fits. 



between the cutoff parameters and the equation of state using simulations. These relations 
turn out to be unaffected by systematics like changes in cosmology (for a fixed equation of 
state) @ and can thus be used to measure the equation of state of the IGM using the cutoff in 
the observed 6(iV)-distribution. 

The amplitudes of the power-law fits to the cutoff and the equation of state are plotted 
against each other in the left panel of Fig. |j. The error bars, which indicate the dispersion 
in the cutoff of sets of 300 lines (typical for z — 3), are small compared to the differences 
between the models. This means that measuring the cutoff in a single QSO spectrum can 
provide significant constraints on theoretical models (at z — 3, physically reasonable ranges for 
the parameters of the equation of state are 10 3 K < T < 10 4 5 K and 1.2 < 7 < 1.7 || §]). 
The slope of the cutoff, Y — 1, is plotted against 7 in the right panel of Fig. [|. The dispersion 
in the slope of the cutoff for a fixed equation of state is comparable to the difference between 
the models. The weak dependence of Y on 7 and the large spread in the measured Y will make 
it difficult to put tight constraints on the slope of the equation of state. 
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